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A recently introduced percolative theory of unipolar organic magnetoresistance is generalized by treating the 
hyperfine interaction semiclassically for an arbitrary hopping rate. Compact analytic results for the magnetore- 
sistance are achievable when carrier hopping occurs much more frequently than the hyperfine field precession 
period. In other regimes, the magnetoresistance can be straightforwardly evaluated numerically. Slow and 
fast hopping magnetoresistance are found to be uniquely characterized by their lineshapes. We find that the 
threshold hopping distance is analogous a phenomenological two- site model's branching parameter, and that 
the distinction between slow and fast hopping is contingent on the threshold hopping distance. 



I. INTRODUCTION 

The prospect of spintronic ff^l [ n organic materialdSEO has 
generated much interest in recent years. Spin-orbit coupling, 
a bane to long spin lifetimes, can often be much weaker 
in organic than in typical inorganic semiconductors used for 
electronics, due to small atomic numbers in the organic con- 
stituents. Affordable manufacturing and chemical tunabil- 
ity also add to the appeal of studying spins in organic sys- 
tems. In contrast to inorganic semiconductors, organic semi- 
conductors are typically disordered and therefore their trans- 
port mobilities are much smaller. Despite this apparent draw- 
back, organic semiconductors are currently used in a variety 
of electronic devices; 7 understanding the behavior of spins 
in these systems offers the possibility of adding magnetic 
functionality to these and future devices. Since spin trans- 
port properties are intertwined with the electrical transport 
properties, 8 features distinct from spin transport through inor- 
ganic semiconductor^^] are expected in organics due to their 
very different electronic transport properties. 

In parallel with these developments, researcheri^L^ have 
studied a magnetic field effect in a diverse array of organic 
materials, the so-called organic magnetoresistance (OMAR). 
Ref. 18 reviews this effect. It exists in nonmagnetic mate- 
rials and is characterized by magnetoresistances of 10-20% 
at fields as small as 10 mT. These properties suggest appli- 
cations in magnetic sensors and organic displays!^ OMAR 
has resisted explanation by typical magnetoresistive effects 
such as Lorentz force deflection, wave function shrinkage, and 
weak (anti-)localization. 21 Several recent models of OMAR 
have been proposed. Most involve spin-dependent processes 
originating from hyperfine interactions and can be classified 
as either unipolar^ or bipola r 12 ! 13 ! depending on whether the 
OMAR mechanism relies on a single carrier type (e-e or h-h 
interactions) or two carriers types (e-h interactions). 

Due to the ubiquitous disorder in organics, methods to cal- 
culate transport properties in inorganic semiconductors fail for 
their organic counterparts. One common technique to eval- 
uate transport in organic materials is percolation 

theoryPSl 

Recently a modeP^has been developed of OMAR which ex- 
plicitly takes into account hopping transport for a single car- 
rier in a disordered material using the theoretical description 
of percolation theory. Here we extend the work of Ref. [24] 
by deriving, from a semiclassical theory of the hyperfine in- 



teraction, similar results that apply within any hopping rate 
regime. 

The percolation model proposed in Ref. 24 and further 
developed here reduces the complex phenomena of spin- 
dependent hopping to a tractable problem of r-percolation 
with an effective density of hopping-accessible sites that is 
modulated by magnetic fields through singlet-triplet transi- 
tions. We focus on unipolar charge transport since several 
analytic results are obtainable; near the end we assess how 
the general features and insights of this model may shed light 
on bipolar magnetoresistance mechanisms. Inclusion of ener- 
getic disorder precludes simple results and is not treated here, 
although similar MR results and trends should be expected. 

OMAR in unipolar transport was studied theoretically by 
Bobbert et al in Ref. [16] and then further developed in Ref. 
[25] within a so-called "two-site" model. In the two-site model 
the resistance is determined by a single "bottleneck" pair of 
sites and a phenomenological branching parameter which al- 
lows carriers to circumvent the bottleneck if the bottleneck 
resistance becomes too large. An additional feature is that the 
two- site model requires a very large electric field to force hop- 
ping in a single direction. Our analysis naturally accounts for 
bottleneck avoidance within percolation theory (no branch- 
ing parameter needs to be introduced) and large electric fields 
are unnecessary. 26 More recently, the same researchers have 
reexamined their two-site model numerically by solving the 
stochastic Liouville equationPS Our model is in qualitative 
agreement with the two- site model on several points as we 
discuss throughout this article. 

Our article is organized as follows: in Section |TT] our the- 
ory is presented; we describe how transport is changed by 
processes that change the relative spin orientation of polaron 
pairs. Section [TTT| sh ows how the MR is calculated from our 
theory. SectionjlV] identifies the hyperfine interaction as the 
MR mechanism and treats it semiclassically. In Section [V| the 
special case of fast hopping is examined because analytic re- 
sults can be obtained. In Section 
arbitrary hopping rates. Section [V 



VI} MR is investigated for 
I examines how our theory 



may pertain to bipolar organic devices. 



II. MODEL 

A spatially disordered organic system can be modeled as 
a network of random resistors as described by Miller and 
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Abrahams for doped inorganic semiconductors. The inter- 
site hopping resistance between two sites, i and j, is given by 
Rij = Roe 2ri j/ £ where r^- is the inter- site separation and t is the 
localization length of a carrier at each site. For simplicity, the 
localization length is taken to be uniform throughout the site 
array. Percolation theory offers a method to calculate the bulk 
resistance in such a random resistor network P 2 ! 29 ! 30 ! A critical 
resistance (distance) R c (r c ), which is the smallest resistance 
(or equivalently the smallest separation) that still allows for 
an infinitely large network of bonds, sets the bulk resistance. 
This percolation length is governed by a bonding criterion: 

B C = A% f C r 2 Ndr, (1) 
Jo 

where N is the density of sites and B c « 2.7 in three dimen- 
sions is a number that determines the average number of bonds 
each site in the percolating network must connect withp^This 
"r-percolative" transport model is valid when the inter- site 
separation is large and temperatures are high, and has been 
observed in organic semiconductors We do not treat here 
smaller inter- site separations and lower temperatures where 
energy disorder is vital to understanding transport. In princi- 
ple, the theory here can be generalized to treat energy disorder. 



triplet polaron pair singlet polaron pair 




FIG. 1. Spin blocking in transport. Top: the initial situation for hop- 
ping a triplet polaron pair (left) and a singlet polaron pair (right). 
Bottom: on the left the spin-blocked carrier has made the more dif- 
ficult further hop to an unoccupied site. On the right, the carrier 
successfully hops to the occupied nearest neighbor since the pair has 
a singlet spinor. The carrier concentration is dilute enough such that 
sites near the polaron pair are unoccupied. 

Figure [T] displays how the Pauli exclusion principle affects 
spin transport in hopping conduction; double occupation at a 
single site is forbidden if their spins form a triplet state (T), but 
permissible if in the singlet state (S)!^^ Following many of 
the earlier works on OMAR,^ we use an alternate terminol- 
ogy from Ref. 24 , and describe carriers in conjunction with 
their localizing sites as polaron quasiparticles. An arbitrary 
spin localized at a site (a polaron) is unable to hop to another 
polaron's site if the polaron pair's (PP) spin state is T, but can 
hop to site if the spin pair forms an S state just as it would to 
an unoccupied site, as schematically shown in Figure [T] The 
respective concentrations of these three types of sites are Nt, 
Ns, and No- Furthermore, spin statistics dictates that Ns = \N\ 



and Nt = ^N\ with N\ being the concentration of singly oc- 
cupied sites. The concentration of carriers is small enough 
such that a polaron encountering a bipolaron is extremely un- 
likely; additionally, PPs are sufficiently separated by unoccu- 
pied sites. 

In physical systems, double occupation of a site costs a 
Coulomb interaction energy U. Within the r-percolation pic- 
ture, U > inhibits double occupation and reduces any spin- 
dependent magnetoresistive effects. However, in more real- 
istic systems when energetic site disorder is larger or on the 
order of U, the effect of U is not as straightforward. In such 
a case, hopping does not occur between sites with identical 
energies but between sites with energy difference U. Sur- 
prisingly, theory for MR in which there is also energy dis- 
order entails larger MR effects for positive non-zero U than 
for U = 0. 34 Since we only consider positional disorder we 
assume U = to avoid an unphysical inhibition of double oc- 
cupation. 




FIG. 2. Spin-blocking in percolative transport. A carrier spin starts at 
the bottom and begins nearest neighbor hops to the top (along solid 
blue line). Hopping-accessible sites with density N' e ^ are denoted by 
black solid circles. Left: Red sites are completely inaccessible due 
to spin-blocking when no spin transitions take place (as described 
in the main text). The resistance is determined by the critical hop 
of r c . Right: when spin transitions occur, some red sites become 
accessible making the total density of hopping-accessible sites N e ff. 
The result of spin transitions is that the average inter- site separation 
and r c decreases which modifies the hopping path. 

Since bipolaron formation is forbidden in the T states, the 
concentration of sites is effectively reduced to N' e ^ —N — Nt 
since only these sites are accessible to a hopping polaron. In 
the absence of T-S transitions (or spin flips in the language 
of Ref. l24k we would then rewrite the bonding criterion of 
Eq. |T|) with N — )> N' e ^. A PP's hopping dynamics is thus 
strongly dependent on T-S transitions since bipolaron forma- 
tion is completely blocked for T-states and allowable for S- 
states. The effective reduction of site density entails that in 
general longer hops need to be achieved and an overall in- 
crease in resistance is expected as shown in Figure [2] (a). If the 
total concentration of singly occupied sites is fixed at N\ , the 
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probability of a successful hop to an occupied site (hopping 
to an occupied site happens with probability N\/N) is 1/4, in- 
dependent of spin effects. So, one-quarter the time the hop 
will be successful and the total density of sites for which suc- 
cessful hops take place is No +Afc. So as before the density of 
unrestricted hopping sites is N' e ^. We must now account for 
the situation that occurs the other three-quarters time in which 
the hopping attempt to a singly occupied site is foiled due to 
occupation by a T- forming spin, which occurs at Nj sites. 

We introduce the possibility that the spin-blocked path can 
be opened by a process that alters the PP's spin state; namely 
transitions from T to S. The probability for the blockade to be 
lifted by the time the next hopping attempt takes place, Xh, is 
Pt^s- We thus modify the effective density of T sites to be 
[1 — pt^s]Nt- The bonding criterion becomes 



B c = 4n J\ 2 N eff dr, 

where further r dependence lies in 

N e ff =N — N T J r Pt^sNt 



(2) 



(3) 



through the hopping time T/^. Our model displays an interplay 
for a PP of two events: waiting for the transition to S to hop 
to the nearest site versus disassociation by hopping to a site 
further away. 



III. CALCULATING THE MAGNETORESISTANCE FOR 
SPATIALLY DISORDERED ORGANIC SYSTEMS 

Eq. ^ is the starting point for deriving the magnetoresis- 
tance. As discussed, the effective density to be used is N e ff, 
which yields 



B c = ^fyi (N-N T ) + 4nfN T [ y y z p T ^sdy, (4) 
3 Jo 

where y c = r c ji is the dimensionless critical length which dic- 
tates the threshold resistance R c = Roe 2yc \ %h = v^ l e 2y is the 
hopping time. A quantity y CQ = (3B c /4n£ 3 N) 1 / 3 is defined as 
the critical inter-site spacing in the absence of all spin effects. 
In general, y c cannot be isolated in Eq. ^ and the resultant 
MR can only be obtained numerically unless the system is in 
the dilute carrier regime (N\ <C N) (which is what is assumed 
throughout this article). 

To solve for the MR we first need to find the critical length 

y c '. 



n£ 3 y 3 c (N-N T )=B C - 4k£ 3 N t y 2 p T ^sdy. 

Jo 



(5) 



3B C 3N T fyc 2 

4%£ 3 {N-N T ) ~ N-N T 



■r 

Jo 



-A 



3N T 



y pr^sdy 

rye 

/ y 2 Pr^sdy. (6) 
Jo 

where y cx = y CQ (l — Nt /N)~ 1 / 3 is the renormalized critical 
inter-site spacing. y c is near y Cl since the singly occupied site 



N-N T 



concentration is small. So on the right hand side we can ap- 
proximate y c ~ y Cl and N — Nt ~ N leaving us with 

/ i 3N T ryci 2 y/3 



^ y 2 Cl N 



1 N T pi 2 



f l y 2 Pr^sdy. 
Jo 



(7) 



We see that by incorporating T^S transitions, the critical 
length decreases from the length where triplet sites are com- 
pletely excluded. R c is then 



2 Nr r^ c \ 1 i 

R c =R e 2yc =R e 2yc ie *i 

o 2 N r r yc i o 

KRoe 2yc l{l __T y 2 pT ^ sdy) . 
3^ci ^ J® 



(8) 



the MR is 



MR 



_R c (B)-R c (0) 



R c (0) 
2 N T 
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T^s(0)-p T ^s(B)}dy (9) 



which was first found in Ref. 24 The problem now re- 
duces to finding the probabilities for singlets at the next hop 
that were initiated in the singlet state. However not all hops 
happen exactly at %h but over an exponential distribution of 
hopping times to account for this, we write pr^s(B) — 
^ Jq pr^siB.^e'^^dt where p T ^ s (B,t) is the density ma- 
trix element signifying the occupation probability of the sin- 
glet state. This quantity pr^s can be related to an easier to 
calculate quantity ps^s which is the population fraction in 
the singlet state that were initially in the singlet state. Their 
relation is |(1 — ps^s(t)). 37 In summary, we find 



2 1 

^2" 



MR = - — 2~T| X 



(10) 



ry c , r°° 1 

y 2 ~ [Ps^s(B) - Ps^s(0)]e~ t/Th dtdy, 
Jo Jo Xh 

where T| =Nt/N. Due to its frequent use, ps^s W1 U be now 
denoted by the simpler p^. To calculate the magnetoresis- 
tance, the singlet population remaining after a time t must 
be determined given various spin interactions. The interac- 
tions considered in the follow sections are the Zeeman and 
hyperfine interactions. Spin- spin interactions (exchange and 
dipolar) are considered elsewhere!^ 



IV. SEMICLASSICAL MODEL WITH NUCLEAR 
MOMENT AT BOTH SITES 



Now the mechanisms by which T-S transitions take place 
are described. The physical picture is that of a PP composed 
of two spin- \ carriers located at two sites. While the spins are 
localized they evolve coherently under the influence of identi- 
cal applied fields and different hyperfine fields. The semiclas- 
sical approximation entails that hyperfine or nuclear spins are 
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accounted for by classical magnetic fields. The hyperfine field 
at a site is composed of many different nuclei as depicted in 
Figure [3j however since the nuclear spin precession is so slow, 
the total nuclear field is assumed constant throughout the po- 
laron's time of residence at that particular site. Since there is 
no nuclear spin order the orientation and magnitude of the to- 
tal hyperfine field varies from site to site. When one of the po- 
larons hops, its coherent spin evolution ceases as the polaron 
now feels a new local magnetic field (if disassociating out of 
the PP by hopping to an unoccupied site) or exists as a bipo- 
laron (if hopping to the other polaron's site) and is necessarily 
in the singlet state (a large exchange interaction prevents fur- 
ther spin evolution due to the proximity of the two polarons). 
If the hopping is fast, the PP spin has had little time to evolve. 
If the hopping is slow, each spin of the PP can be thought of as 
having precessed many times around its local magnetic field. 

The mathematical details of the semiclassical approxima- 
tion now follow. The Larmor frequency of a polaron spin lo- 
calized at a site due to the nuclear conglomerate spin is the 
constant classical vector 



(11) 



where aj is the hyperfine coupling constant in angular fre- 
quency units between the electron and the j-th nucleus. In 
is constructed from many nuclei each with vector length 
aj + 1) pointing in a random direction. Ij is the spin 

quantum number. The probability distribution for finding a 
specific site among an ensemble of such sites with its total 
end-to-end vector between + dl^ i JSEH 

W(I N ) = (47r4)- 3 / 2 exp [-(4/44)] (12) 




FIG. 3. Semiclassical description of total field. Each site is com- 
posed of some number of spin-moment-carrying nuclei which com- 
bine with the external field GDo to create (0^. Each site has a different 
total field picked from the distribution of Eq. (16) . 



distribution W. In Section |VI| for computational reasons, the 
scheme of Ref . EH is followed by integrating over the total 
field (0%) as opposed to the hyperfine field. To aid this the 
probability distribution is rewritten as 



/ 1 \3/2 / 



1 cog + c4 - 2co cOiv COS 0^v 

4 air 



(16) 

with differential volume element sin 0^ d 0^0)^ d COw d . To 
find the MR, determinations of the singlet density matrix ele- 
ments must be made. 



where 



i). 



(13) 



The effective hyperfine coupling due to all the nuclei at a site 
is a e ff = ly/la^ (although other conventions do exist 40 ). a e ff 
could be different for different types of molecular sites in 
which case it would have to be labeled by a site index (we ne- 
glect this effect here). The effective hyperfine magnetic field 
is B e ff = a e ff/y e with the electronic gyromagnetic ratio being 
y e = 0.176 ns _1 mT _1 . The effective magnetic field is on the 
order of 1 mT or a e ff ~ 0.176 ns _1 for many organic systems 
demonstrating OMAR. 

The total precession rate seen by a carrier at a single site is 
then 



(0^=7^ + 0)0, 



(14) 



where <0o = J e Bz is the applied field. So the PP Hamiltonian 
is 



] N 2 



•s 2 , 



(15) 



where site indices have been incorporated. To account for the 
ensemble of PPs the singlet probability is averaged over the 



V. FAST HOPPING 

Before examining a theory applicable to any hopping time 
and field strength, we first explore the fast hopping regime for 
which simple and tractable analytic results can be obtained. 
The results derived herein allow us to confirm the validity of 
the spin relaxation model proposed in Ref. 24 over the entire 
range of magnetic fields. 

The dynamics for small a^t can be solved analytically 
for arbitrary cOq^ using perturbation theory in the interac- 
tion representation!^ First Eq. (15) is rewritten as = 
Mhf + Mz where 



3^ 



Im - Si +^v 2 S 



and 



hf =I Ni 



) So 



(17) 



(18) 



are the hyperfine and Zeeman Hamiltonians, respectively. In 
the interaction representation, the following operators are de- 
fined: 



p*(0 



P(0 



(19) 



(20) 
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initially it can be shown that p*(0) = p(0). Time depen- 
dent perturbation theory for the density matrix to second order 

gi ves l43J 

P , (0 = p(o) + ^[p(o),^ / ( ( 'p'- 

hfoll M )^h}(t%^h}(t'')W'dt'. (21) 



The new hyperfine Hamiltonian is 

J? h * f (t) = U(Jf hfi +J? h f 2 )U t , 

with 



U 



® t 



(O t 



(22) 



[cos + 2iS\ sin —^-} [cos + 2iS\ sin —^-\ , 

(23) 

because Si commutes with S2. In our singlet/triplet basis 
we write p(0) = Ps, which is the singlet projection operator 
written in Appendix [A] The singlet part of the density matrix 
(S|p*(f)|S) is ps(t). Initialization in the singlet state requires 
that ps(0) = 1. The first order correction vanishes. After av- 
eraging over the hyperfine fields for the ensemble of two car- 
riers, including the second order term yields 



1 

16 



1+2 



sin 2 0)0/72 



(<W/2) 2 j 



(24) 



in agreement with the quantum mechanical calculation ! 42 ! 44 ! 
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FIG. 4. Fast hopping normalized magnetoresistance for y c = 5. Red: 
r = 0.5 x 10 6 ; blue: tq = 1 x 10 6 ; green: tq = 2 x 10 6 . Accom- 
panying points are calculations from general theory of Section |Vl] 
Not shown in plot is that the saturated MR strongly reduces with an 
increased hopping rate [which can be seen from Eq. (27)]. 

To find the MR, we substitute the singlet probability, Eq. 



( [24] ), into Eq. ( [10) . After integrating over the exponential 
distribution of hopping times, 



MR=-\\- r — ; J / 2 y l dy. (25) 

6 3^ Jo 



co 2 + l/x 2 



This results agrees with the spin relaxation result 24 to within 
a numerical factor. The agreement validates the interpretation 



given by the spin relaxation model; the intersystem crossing 
between singlets and triplets can be described by spin relax- 
ation due to the rapidly varying hyperfine interaction (due to 
the fast hopping) in the motional narrowing regime. The spin 
relaxation rate is taken to be 



= a eff\ 



7(co§ + x,- 2 ). 



(26) 



The probability for the spin-flip is 1 — exp(— Th/i s ) ~ 

The integral in Eq. ( [25) can be computed when the hopping 
rate has an exponential dependence on the hopping distance 
\ l = vo exp(— 2y). The result is cumbersome but can be con- 
siderably simplified under the usual assumption that y c ^> 1 
to 



MR=2^(a 2 //;1 



1 



ln(l 



, (27) 



where x c 1 = voexp(— 2y c ). At large fields, the MR saturates 
at MR sat = ^r\(al ff l +a 2 // 2 )x 2 . Figure H shows several in- 
stances of the fast hopping MR. The following general fea- 
tures should also be noted. First, the saturated MR is depen- 
dent on the hyperfine field and the hopping rate; the MR de- 
creases as hopping times shorten, as the hyperfine fields have 
less time to mix the triplet to the singlet state. Second, the 
MR lineshape is independent of the hyperfine field and solely 
dependent on the hopping rate. The width increases with in- 
creases in the hopping rate; larger fields are required to sup- 
press the hyperfine fields as evident from the motional nar- 
rowing spin relaxation formula Eq. ( [26) . Finally, the derived 
MR is always positive; the applied field suppresses T-S mix- 
ing. These features were pointed out in Ref. |24|and were also 
confirmed by numerical simulations solving the stochastic Li- 
ouville equation!^ 



VI. ARBITRARY HOPPING 

The results of the previous section are valid only for short 
hopping times. A different approach must be utilized to evalu- 
ate the MR for long hopping times. Initially the density matrix 
is p(0) = \S) (S\. At some later time, under the evolution of the 
Hamiltonian, J$f, 



p(t) = cxp(-iJ^t/h)p(0)Qxp(iJ^t/h) 
= exp(-ijet/h)\S)(S\exp(iJt?t/h). 



(28) 



We are interested in the S portion of the density matrix ps = 
(S|p(f)|S) which then becomes 

ps = (S\exp(-iJ$?t /h)\S)(S\exp(iJ?t /h)\S) (29) 
= \ (S\exp(-ijet/h)\S)\ 2 . 

An average over the nuclear field distribution is taken to ac- 
count for the ensemble of carriers at sites with differing nu- 
clear configurations. Under certain restrictions for the n uclei, 
the problem can also be solved quantum mechanically!^^ 
However when many nuclei are present at each site, which is 
the case in disordered organic semiconductors, the quantum 
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mechanical calculation is best tackled numerically. As ex- 
pected, it has been shown that the validity of the semiclassical 
approximation improves with the number of nuclei!^ 



Eq. ( 29 ) is solved by noting that our Hamiltonian, Eq. ( 15 ) 
is Zeeman-like so we can use: 



-icn-S 



c c 
- cos - — 2in • S sin - . 
2 2 



(30) 



In Eq. ( |3Q| ), we have the total field at a site (G>#) unit vector 

n x = sin0^cos(|)iv, n y = sin 0^ sin (|)#, n z = cos0^, (31) 

with angles defined in Figure [3] With some labor, it can be 
shown maPS 



Ps(B = 0) 



(38) 



Refs. [35] [44] and |40] provide generalizations for when the 
two sites are of different hyperfine species. For slow hopping 
we find that it is favorable to perform the time integral of Eq. 
( [36] ) first which can be done analytically though we omit it 
here because the expression is cumbersome. The integral over 
u is then conducted numerically. Nevertheless the extreme 
hopping MR converges to the following expressions: 



MR^(vo^O) 



1 

27 3 



(39) 



Ps = Fi(l)Fi(2) + 2F 2 (1)F 2 (2) + F 3 (1)F 3 (2) + -F 4 (1)F 4 (2), 

(32) 

where 



2 C0^ v „ ,. x 1 



= 1 - <sin z (^));, F 2 (i) = -(n z (i)sm(a N t))i (33) 



F 3 (i) = (^(/)sin 2 (^)) ; -, F 4 (i) = 1-F l (i)-F 3 (i) (34) 

where i refers to site one or two and angular brackets refer 
to averaging over total field distribution W. There are three 
unique averages that need to be calculated: 

h(i) = (sin 2 ( C y^)),-, I 2 (i) = ^(y N shx(<o N t))i, 



/ 3 (0 = (^sin 2 (^)), 



(35) 



where we have made the change of variable = cos0^. Eq. 
( [32] ) can be expressed in closed form though we refrain from 
doing so for the sake of brevity. We still need to integrate over 
time and radius: 

MR=-— x\ [p 5 (B)-p 5 (0)] / y 2 -e-'^dydt, 
3^ Jo Jo 

(36) 

which is never a negative value if hyperfine coupling widths 
are identical. In general results are achieved by numerical in- 
tegrals over y and t. We confirm that our general calculation 
agrees with the analytic results of Section [V] (solid symbols 
in Figure]?]). However we find that performing both integrals 
numerically is most practical in the intermediate to fast hop- 
ping cases. We find that in the slow hopping regime, making 
a change of variable u = exp(— 2y) improves the ease of nu- 
merical evaluation. 



A. Saturated Magnetoresistance 

The singlet probabilities simplify at zero field and infinite 
fields. Hence it is instructive to examine the saturated MR. 
The singlet probabilities reduce to the following: 



1 . -^//4 



(37) 



MR^(v ^°o)=0, 



(40) 



which qualitatively agrees with the simpler model of Ref . [24] 
Figures [5] depicts the saturated MR versus the hopping rate. 
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FIG. 5. Saturated magnetoresistance as a function of hopping rate. 
The inset is same as main but focuses on fast hopping rates by plot- 
ting on log-log graph. Red: y c = 5; blue: y c = 7. Solid lines are 
Taylor expansions around zero and infinite hopping rate with 50 
terms. Note that the slow hopping regime extends considerably past 
vo/a ef f = 1. 

The overall shape is similar to Ref. [27] though the decrease 
occurs at larger vo here (to be discussed below). Also in con- 
trast, the shape is slightly non-monotonic before the sharp de- 
cline. The hopping rate dependence highlights three interest- 
ing features: 

Correspondence between critical radius and branching ra- 
tio. The critical radius y c acts analogously to the two-site 
model '^2^ branching ratio b = r a p/r a?e where r a p is the 
rate from occupied site a to occupied site P and r a , e is rate 
from occupied site a to the environment (in essence, avoiding 
the occupied site (3). A high branching ratio entails that the 
polaron spin's only way to move off of a is to hop to (3 (and 
can only do that if they are a singlet pair). Our critical ra- 
dius acts similarly; large y c entails small site density and large 
inter-site spacings. Since sites are so far apart and the hop- 
ping rate decreases exponentially with distance, if the nearest 
site happens to be spin-blocked, the polaron at a will likely 
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wait until the spin configuration is favorable instead of the ex- 
tremely difficult further hop to a next nearest neighbor (anal- 
ogous to the environment of the two-site model). Hence hops 
to P occur more frequently than hops not to P, much like is 
phenomenologically modeled by the branching ratio param- 
eter in the two- site model. In cases of large y c (or b), the 
saturated MR is larger, as the carrier spin must depend solely 
on the T-S transition (and hence the applied field) — there is 
no possibility to avoid the occupied site. 

Transition from fast to slow hopping occurs at an unex- 
pected hopping rate. The terms 'fast' and 'slow' hopping are 
not as simple as to define as vo ^> 1 and vo <C 1, 
respectively. This is because the spatial dependence is also 
important and effectively decreases the hopping rate. Also 
voexp(— 2y c )/a e ff = 1 is not a good measure of the criterion 
because most hops occur across distances less than y c . There- 
fore we predict that the "slow hopping" regime is in fact ap- 
plicable at faster hopping rates (vo) than previously expected. 
This prediction is consistent with experimental observations 
that large OMAR occurs with hopping rates expected to be 
faster than the hyperfine frequency!^ 

Limiting cases of the saturated MR for slow hopping. When 
hopping is slow, the saturated MR is independent of both the 
hyperfine coupling and the hopping rate. This result is sensi- 
ble since after long waiting times, the PP spins have sufficient 
time to fully mix. 



B. Slow Hopping Magnetoresistance Curves 

The magnetoresistance lineshapes are not fundamentally 
different from the discussion in Ref. |24| Figures [6] and [7] show 
MR traces calculated for several hopping rates at the threshold 
radii y c = 5 and y c = 1. 



20 




5 10 15 20 



Bo/B & ff 

FIG. 6. Magnetoresistance at y c = 5. Black: vo/a e ff = 5 x 10 4 ; 
red: vo/a e ff = 1 x 10 4 ; blue: vo/a e ff = 1 x 10 3 ; green: vo/a e ff = 
1 x 10 2 ; orange: vo/a e ff = 1 x 10 1 . The magnetoresistance is an 
even function of Bo/B e ff. 

In the slow hopping regime (as depicted in Figure [5}, the 
MR width is independent of hopping rate in sharp contrast to 
the fast hopping case (see Figure]?]); the width varies linearly 
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25 




FIG. 7. Magnetoresistance at y c = 1. Black: vo/a e ff = 1 x 10 6 ; 
red: vo/a e ff = 1 x 10 5 ; blue: vo/a e ff = 1 x 10 4 ; green: vo/a e ff = 
1 x 10 3 . The magnetoresistance is an even function of Bo/B e ff. 

with the hyperfine coupling strength which also is a different 
behavior than seen in the fast hopping regime. The large MR 
widths (~ 40 mT ^> B e ff) measured by some researche d 15 l 47 l 
suggest that those scenarios were fast hopping where the hop- 
ping rate's role in MR width is indeed present. 

MR is always positive, in disagreement with the ultra-small 
field effect observed in the simulations of Refs. I48landl27l 
At this time the source of the discrepancy between our results 
and their simulations is not known. As discussed throughout 
this article, on other points the two approaches are in qualita- 
tive agreement. Positive MR (ignoring ultra-small field effect) 
has been observed in unipolar diodes. 49 It is noteworthy that 
spin- spin inte ractions also cause an ultra-small field effect to 
occur ! 27 ! 38 ! 49 ! Additionally if nuclear spin moments are con- 
sidered quantum mechanically, an ultra- small- field effect is 
expected as shown by Ref. [49] for a single nucleus. We do not 
expect a small number of nuclei per site in the organic systems 
considered here so the semiclassical approximation is valid!^ 

Our result suggests that organic materials with large y c 
(small localization length or small site concentration) yield the 
largest MR. Though increasing bias voltage tends to increase 
the localization length, 50 and therefore decrease MR in our 
theory, experimental observations 18 of the bias dependence 
are unclear since majority and minority current injection rates 
possess a bias dependence for a bipolar organic device. 



VII. APPLICABILITY TO BIPOLAR SYSTEMS 

Our theory predicts solely positive MR. This is simple to 
understand since applying a field suppresses hyperfine in- 
duced T —> S transitions which causes a carrier to either wait 
for the transition or make a slower hop. This slowing down of 
carrier hopping leads to the increase in resistance. 

However the majority of experiments have observed nega- 
tive MR. 14 18 This is due to the presence of two types of car- 
riers (bipolar system). While the details of our theory do not 
apply in bipolar systems, we can still see qualitatively why 
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negative MR might be dominant in a simplified model of a 
bipolar organic device. 

Two oppositely charged polarons (an exciton) at the same 
site do not contribute to the current since they will either re- 
combine (luminesce) if a singlet or remain as an exciton if 
a triplet due to the large attractive Coulomb interaction (ex- 
citon disassociation is ignored for simplicity). This is very 
different than the unipolar case we considered where easy for- 
mation of bipolarons encourages current flow. If exciton for- 
mation varies between singlet and triplet e-h pairs, a similar 
spin-blocking mechanism emerges; in dramatic contrast to the 
unipolar scenario, this time more spin mixing leads to more 
exciton formation which inhibits current. An applied field 
suppresses spin mixing (again only considering the hyperfine 
interaction) which leads to less exciton formation, more cur- 
rent, and therefore negative MR. Developing a quantitative 
theoretical framework for bipolar OMAR based on percola- 
tion theory is a challenge for future investigations. 



VIII. CONCLUSION 

Testing our theory quantitatively is most tractable for high 
temperatures and a low density of molecular sites. Controlling 
the density of sites is predicted to change the hopping rate and 
the resulting MR. Such a manipulation of site densities has 
been successfully employed in the past in TNF films in which 
conduction via r-percolation was measured through time-of- 
flight experiments. 31 32 In these experiments the molecular 
density of TNF was carefully controlled by dispersing TNF 
in an inert polyester host. 

The theory presented here has implications for MR effects 
in amorphous semiconductors, El colloidal quantum dotsp2 
spin diffusion in organic full spin valves, 53 and MR effects in 
organic semi- spin valves where fringe fields from a magnetic 
film create a unique MR curve. El 
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Appendix A: Spin operators 

We write all matrices in the singlet/triplet basis. The spin 
ladder operators are 




and S t = . The other spin operators are = S * ^J~* 



and 



S z — - 



/o 1 
1000. 

1 I ' 2 2 
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/ -1 
-10 
10 
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The singlet projection operator is 
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Appendix B: The Ii,h, and I3 integrals 

There are three unique integrals to calculate. First, 
1 



h = 



4% l / 2 a 3 N 



1 

e nx 



(Bl) 



f°° 9 "4T 9/ 

/ toid(d N e a N sin z/ 
Jo 



1 gwyyjv 



(-2~) J^Ne 
or in dimensionless units xn = co/v/a/v 



4jiV2 



(B2) 



rx 2 N dx N e-* x * sin 2 (*^) C dy N e^^. 
Jo 2 J -1 



which is 



11 x 

I 1 = -+[- cos(/¥c) - — sin(Ax)] e~ 



2 L 2 v ; Sh 
where h = (tio/a e ff and x = a e fft. Also 
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The last integral is 
1 
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with the result 
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h=h 



: D(V2h)- 



-2h 2 



\6V2h 3 



v / S(Erf( 



2V2 



iy/2h)-Erf( 
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where D(z) = e z J^e x dx is Dawson's integral and Erf(z) 
^ Jq e~ x dxis the error function. 
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